The first-principles study of thermodynamical properties of random magnetic 

overlayers on fcc-Cu(OOl) substrate 
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We present the theoretical study of thermodynamical properties of fcc-Cu(OOl) substrate covered 
by iron-cobalt monolayer as well as by incomplete iron layer. The effective two-dimensional Heisen- 
berg Hamiltonian is constructed from first principles and properties of exchange interactions are 
investigated. The Curie temperatures are estimated using the Monte-Carlo (MC) simulations and 
compared with a simplified approach using the random-phase approximation (RPA) in connection 
with the virtual-crystal approach (VCA) to treat randomness in exchange integrals. Calculations 
indicate a weak maximum of the Curie temperature as a function of composition of the iron-cobalt 
overlayer. While a good quantitative agreement between RPA- VCA and MC was found for iron- 
cobalt monolayer, the RPA- VCA approach fails quantitatively for low coverage due to the magnetic 
percolation effect. We also present the study of the effect of alloy disorder on the shape of magnon 
spectra of random overlayers. 

PACS numbers: 75.70.Ak, 75.50.Bb, 75.30.Et, 75.40.-s, 05.10.Ln 
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I. INTRODUCTION 

The Curie temperature is one of the most important 
characteristics of ferromagnetic materials. Its parameter- 
free determination for bulk ferromagnets and their alloys 
has progressed in the last decade. A reasonable agree- 
ment between calculated and experimental Curie tem- 
peratures was found for transition metal ferromagnets, 
some ordered and disordered transition metal alloys (e.g., 
Ni-based fcc-alloys), some f-metals (e.g., hcp-Gd, bcc- 
Eu), diluted magnetic semiconductors (e.g., GaMnAs) or 
Heusler alloys (see, e.g., a recent review)^. This progress 
was due to a combination of the first-principles determi- 
nation of parameters of a (classical) Heisenberg Hamil- 
tonian and its study using sophisticated statistical meth- 
ods like, e.g., the random-phase approximation (RPA) 
or Monte-Carlo (MC) simulations. It should be noted, 
however, that there are magnetic systems for which the 
above approach has certain limitations or even fails. The 
present approach assumes the existence of robust mo- 
ments and dominating effect of pair exchange interac- 
tions. Systems with induced moments, (e.g., FeRh), mag- 
nets with more complex than pair interactions, etc., are 
just few such cases, where the above approach is not 
successful. It is known that magnetic moments at sys- 
tem surfaces are enhanced due to the reduced number of 
nearest-neighbors, which is favorable for the validity of 
the Heisenberg model. 

Determination of Curie temperatures Tq of low- 
dimensional systems such as ultrathin films or even ran- 
dom magnetic monolayers deposited on non-magnetic 
substrates was studied very rarely in the pasl^^ de- 
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spite their importance. Regardless of considerable efforts 
in the past decade, the parameter-free determination of 
the Curie temperature of low-dimensional systems in the 
framework of itinerant magnetism remains a challenge 
for the theory. There are several reasons for this: (i) 
Already determination of the electronic structure of low- 
dimensional systems, in particular on the first-principles 
level, is much more demanding as compared to bulk sys- 
tems, (ii) the presence of randomness in the system is an- 
other complication, (iii) the statistical treatment of the 
two-dimensional systems is much more delicate problem 
as compared to bulk systems because in two dimensions 
the interactions decay with the distance more slowly than 
in the bulk, and (iv) the presence of relativistic effects has 
to be taken into account for monolayers (the Mermin- 
Wagner theorem). 

This paper is a natural extension of our previ- 
ous study^ in which the Curie temperature of two- 
dimensional systems was formulated for ideal, non- 
random systems based on exchange integrals deter- 
mined from first-principles (specifically for Fe- and Co- 
overlayers on the fcc-Cu(OOT) substrate). Here we wish to 
extend this study to random overlayers. As case studies 
we consider (i) The random (Fe,Co)-overlayer deposited 
on the fcc(OOl) face of Cu, and (ii) The fcc-Cu(OOl) sub- 
strate with an incomplete coverage by iron atoms. One of 
aims of this study is a development of a reliable scheme 
for determination of Tq which employs the MC method. 
This approach is accompanied by a simplified approach 
based on the RPA method in which randomness in ex- 
change integrals is treated approximately in the virtual 
crystal approximation (VCA) and its reliability and limi- 
tations are tested. We also present the study of the effect 
of alloy disorder on the shape of magnon spectra of ran- 
dom (Fe,Co)-overlayers obtained by spin dynamics. 

There are no experimental data for the present sys- 
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tern to compare with although FeCo-overlayers on var- 
ious non-magnetic substrates (Cu,Pd,Rh) were studied 
often with the aim to estimate the magnetic anisotropy 
energies 4-6 . Also, there are no doubts on the importance 
of a parameter- free approach to estimate Tq of imperfect 
ultrathin magnetic overlayers (see, e.g., a recent review^). 



II. FORMALISM 
A. Electronic structure and Heisenberg model 

The electronic structure of the system was determined 
in the framework of the Green function implementation of 
the scalar-relativistic tight-binding linear muffin-tin or- 
bital method (TB-LMTO) in which the effect of the semi- 
infinite substrate was included properly in the framework 
of the surface Green function (SGF) approach while the 
disorder in the overlayer was treated in terms of the co- 
herent potential approximation (CPA). The case of par- 
tial coverage of fcc-Cu(OOl) by iron atoms is simulated as 
a random overlayer consisting of iron atoms and vacancies 
described by empty spheres. The vacuum above the over- 
layer was simulated by empty spheres (ES). Electronic 
relaxations were allowed in four empty spheres adjoining 
the overlayer, the overlayer itself, and in five adjoining 
Cu substrate layers. This finite set of layers was sand- 
wiched selfconsistently between a frozen semi-infinite fec- 
Co(OOl) and the semi-infinite vacuum including the effect 
of the dipole surface barrier. Possible small layer relax- 
ations between overlayer and substrate were neglected 
although the present approach allows to include them 8 . 
We refer the reader for more details to Refs. loU 1 OL 

An important advantage of the TB-LMTO-SGF ap- 
proach is a possibility to estimate exchange interac- 
tions between magnetic atoms in the overlayer by a 
straightforward generalization of the well-approved bulk 
conceptiii 1 -. The exchange integrals jff* between 
sites i,j occupied by atoms Q and Q' (Q,Q'=Fe,Co or 
Q, Q'=Fe,vacancy) in the magnetic overlayer may be ex- 
pressed as follows 1 - 
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Here, the trace extends over s,p, d— basis set, the quanti- 
ties df are proportional to the calculated exchange split- 
ting, and the (auxiliary) Green function describes 
the propagation of electrons of a given spin (a =f, I) be- 
tween sites i, j in a random overlayer (the bar denotes the 
CPA configurational averaging) . The integration path C 
in the complex plane starts below the bottom of the va- 
lence band and ends at the Fermi energy. It should be 
noted that both the direct propagation of electrons in 
the random magnetic overlayer and the indirect one in 
the semi-infinite Cu-substrate are included in Eq. (fT]) on 
an equal footing. Finally, the disorder-induced vertex 
corrections due to the correlated motion of two electrons 



in a random overlayer can be neglected in a reasonable 
approximation due to the vertex-cancelation theorem 1 ^. 

Once the exchange interactions were known, we con- 
structed a two-dimensional (2D) random classical Heisen- 
berg Hamiltonian to describe the magnetic behavior of 
the random (Fe,X)-overlayer (X=Co, vacancy) on a non- 
magnetic fcc-Cu(OOI) substrate 

Q,Q'=Fe,Co ijtj 

+ £ A«5>?) 2 (ef) 2 . (2) 

Q=Fc,Co i 

In Eq. @, e.i denotes a unit vector with the direction 
of the local magnetic moment at the site i and r]f is the 
occupation index which equals 1 if the site i is occupied 
by the atom Q and zero otherwise. The second term is an 
uniaxial anisotropy of strength A® with an easy axis out 
of plane. By construction, the value of the corresponding 

magnetic moment is included in the definition of , 
and positive (negative) values denote FM (AFM) cou- 
plings. Because exchange integrals for two-dimensional 
case decay with the distance d between sites i and j more 
slowly as compared to the bulk, a large number of shells 
is needed to obtain well converged results, in particular 
for non-random cases of pure Fe- and Co-overlayers. Up 
to 90 shells were included in calculations of Tq in the 
framework of the RPA and MC methods. The effect of 
small induced moments on substrate atoms was neglected 
(moments are of order few hundredths of /xb). 



B. Statistical treatment 

1. Random Phase Approximation 

The expression for Tq in the RPA is a generalization 
of its bulk counterpart 1 to the case of random magnetic 
overlayers 2 : (i) A vanishing Tq is obtained in agree- 
ment with the Mermin- Wagner theorem for vanishing 
anisotropy energy ; (ii) The anisotropy energy is taken 
here as an adjustable parameter which we have identified 
with the dipolar energy and used the same values as in 
Ref. l2|. This is not a serious problem as Tq has only a 
weak logarithmic dependence on A Q , cf. Ref. [H; and 
(iii) We have averaged three exchange integrals 
(Q, Q'=Fe,Co or Q, Q'=Fe,vacancy) and introduced the 
effective non-random exchange integrals Jff defined as 
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which depend on the actual composition of random alloys 
Fe x Xi^ x (X=Co, vacancy). This is the virtual-crystal 
approximation (VCA) and we have tested its applicabil- 
ity by performing MC simulations in which this approx- 
imation is not used (cf. Sec. MI B[) . Using the VCA, the 
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Curie temperature is 
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where J cS (q\\) is the lattice Fourier transform of the ef- 
fective exchange integrals Jff . 



2. Monte Carlo simulation 

In present work, we use the UppASD package^, 
developed at the Uppsala University for the study of 
magnetic materials. This package also contains MC 
code with the Metropolis algorithm which we utilized. 
The sampling in the Metropolis algorithm is controlled 
by the transition rate w(i, /) from an initial state i 
to a final state / which depends on energies 2% and 
Ef of an initial and a final state, respectively. The 
energies Ei and Ef were calculated using the Heisenberg 
Hamiltonian ©. 

The number of trial moves in every Monte Carlo step 
(MCS) corresponds to the number of spins in the over- 
layer. The substrate is treated as rigid, adatoms do not 
diffuse, but rotation of spins of all adatoms is allowed. 
We start the simulation with a system of randomly dis- 
tributed spin vectors in the overlayer. In the first phase 
of simulation, we only equilibrate the spin system; we 
use typically 10,000 MCS. After that, we perform the 
measurement of observable quantities; we measure over 
50,000 MCS. We average calculated quantities over sev- 
eral independent runs. In the case of a random alloy, we 
employ several (typically 20) different random distribu- 
tions of Fe and Co adatoms or Fe adatoms and vacancies 
in the overlayer. 

One can use several methods for an estimation of 
the Curie temperature. We illustrate these possibili- 
ties in the case of a random (Fe,Co)-overlayer with the 
concentration x = 0.5. One option is to employ the 
fourth-order size-dependent Binder cumulantM JJ L = 

1 — (m*) I 3 (ml) 2 . Here, m z is the z-component of 

the magnetization and L is the system size. This ap- 
proach turned out to be quite useful in the past, though 
some limitations exists. The shape of temperature de- 
pendence of the magnetization is influenced by the size 
L. It leads often to the situation that cumulants for dif- 
ferent system sizes cross in one point corresponding to 
the Curie temperature^. We calculated cumulants in 
our case, see Fig. [TJ Unfortunately, we have found that 
the identification of a crossing point is not so clear as in 
the case of three dimensional systems. 

Other methods are based on finding of a singularity 
point in the susceptibility \ or in the specific heat C. 
In the simulations, one can directly obtain only size- 
dependent quantities xl or Cl- One locates the tem- 
perature T£ or Tq of a local maximum of the suscepti- 
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FIG. 1: (Color online) Temperature dependence of Binder 
cumulants Ul for different system sizes in the case of fee 
(Fe,Co)/Cu(001) random overlayer with the concentration 
x = 0.5. 



bility or specific heat for the size L and then performs 
extrapolation using corresponding scaling relation. 

In this paper, we used mainly susceptibility (see Sees. 
IIIIBI and MI Cp . After location a temperature of a 
local maximum of the susceptibility for the size L, we 
utilized the scaling relation 



(5) 



where T x is an estimate of Curie temperature for the 
infinite system. More specifically, the critical exponent 
v = I for the two dimensional Heisenberg model with 
uniaxial anisotropy is known 16 . The system falls into 
the same universality class as the two-dimensional Ising 
model, for which all the critical exponents have analyti- 
cally known values. 

For the comparison with the method of cumulants, we 
evaluated Tq by the calculation of the size-dependent sus- 
ceptibilities xl for fee (Fe,Co)/Cu(001) random overlayer 
with the concentration x — 0.5 in Fig. [5J We carried 
out simulations for several system sizes (typically rang- 
ing from L = 16 to L = 128 and then we estimated Tq as 
an extrapolation of Tq using linear regression for infinite 
L. 

Similarly we calculate the size-dependent specific heat 
(figure is not presented). 



C. Atomistic spin dynamics 

Using the generalized Hamiltonian H , Eq. @, as a 
starting point, the temporal evolution of the atomic mo- 
ments, mi, where mi = |mi|ei and \rrii\ is the amplitude 
of the magnetic moment, at finite temperature is gov- 
erned by Langevin dynamics through coupled stochastic 
differential equations of the Landau-Lifshitz form, 
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FIG. 2: (Color online) Temperature dependence of size- 
dependent susceptibilities \l for fee (Fe,Co)/Cu(001) random 
overlayer with the concentration x — 0.5. 
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dt (1 + a 2 ) 
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m, x [Bi+bi(t)] 

m, x {m. t x [B; + bi (*)]}, (6) 



m|(l + a 2 ) 



where 7 is the electron gyromagnetic ratio, and is 
the effective field produced by all other moments. Ef- 
fects of temperature are included via adding a stochastic 
Gaussian-distributed magnetic field b,(£), the strength 
is controlled by the Gilbert damping parameter a. In 
principle, the Gilbert damping parameter could be cal- 
culated by means of electronic structure methods, but 
here it is considered as a parameter. More information 
about the numerical details and integrating scheme for 
solving Eq. (j6|) can be found in Ref. [T7|. 

The principal advantage of combining first-principles 
calculations with the atomistic spin dynamics (ASD) ap- 
proach is that it allows to address the dynamical prop- 
erties of spin systems at finite temperatures^— . Two 
important quantities we focus on are the space- and time- 
displaced correlation function, 



C£(t) = «(t)m?(0)) - «(i))<m£(0)), (7) 

where the angular brackets denote an ensemble average 
and v is the cartesian component. Its Fourier transform 
is the dynamical structure factor 



S u {q ]U ) = 
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e**C%(t)dt, (8) 



where q and w are the momentum and energy transfer, 
respectively. S*(q, uS) is the quantity probed in neutron 
scattering experiments of bulk systems, and can analo- 
gously be applied to spin-polarized electron energy loss 



spectroscopy (SPEELS) measurements. By plotting the 
peak positions of the structure factor along particular di- 
rections in reciprocal space, the magnon dispersions may 
be obtained— ~—. 



III. RESULTS 

A. Exchange interactions 

In this section we wish to illustrate some general fea- 
tures of exchange interactions for 2D systems. We shall 
start with the study of the dependence of exchange in- 
tegrals jJ°' Fc on the distance between Fe atoms in the 
fcc(OOl) monolayer along the direction [1,1] in various 
cases, namely (i) The isolated (unsupported) Fe-layer, 
(ii) The Fe-overlayer on fcc-Cu(OOl), and (iii) The ran- 
dom (Feo.5,Coo.5)-overlayer on fcc-Cu(OOl). In this way 
we can study the effect of indirect interactions of Fe- 
atoms via the substrate (missing for unsupported layer) 
as well as the effect of disorder. The case of unsupported 
Fe-layer was studied using the same model in which, how- 
ever, the iron overlayer was separated from the substrate 
by eight layers of empty spheres from the fcc-Cu(OOl) 
substrate. It is well-known that the exchange integrals 
in bulk ferromagnets decay with distance d as d~ 3 . It 
is seen from Fig. [3h. that for the unsupported layer the 
decay with distance is much slower, namely, proportional 
to d~ 2 . 

The effect of the substrate is striking as illustrated in 
Fig. [3p for Fe-overlayer on fcc-Cu(OOl). The decay of 
exchange interactions with distance d is approximately 
d~ 2 5 as a result of interactions via the substrate. Also 
illustrated in Fig. is the additional exponential decay 
of exchange interactions due to the alloy disorder for the 
case of (Feo.5Coo.5) overlayer on the substrate. 

The exchange integrals for pure Fe and Co overlay- 
ers on fcc-Cu(OOl) substrate are shown in Fig. We 
observe generally larger values of exchange integrals for 
the iron overlayer as compared to a cobalt one indicat- 
ing a higher Curie temperature in the former case. Also 
shown are configurationally averaged integrals , 
Q, Q'=Fe,Co for (Feo.5,Coo.5) overlayer (for VCA values, 
see Eq. ([3])) used in an approximate treatment of the 
Curie temperature (see below). It is interesting to note 
that these effective integrals are similar to the species 
resolved integrals shown in Fig. |4p. 

Species-resolved exchange integrals for equiconcetra- 
tion coverage of fcc-Cu(OOl) substrate by iron and cobalt 
atoms are shown in Fig. |4Jd. The remarkable feature is a 
relative similarity of all three kinds of exchange integrals 
which is an indication of the validity of a simplified VCA 
treatment as discussed below. 

Finally, in Fig. [5] we show the dependence of exchange 
integrals J?' Fe on the coverage x of the fcc-Cu(OOl) sur- 
face by iron atoms. The exchange integrals increase with 
decreasing coverage. The reason is their impurity char- 
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FIG. 3: The exchange integrals as a function of the distance 
along the direction [11] in the magnetic layer: (a) The un- 
supported Fe- layer, and (ii) The case of Fe- and Feo.sCoo.s 
overlayers on the fee Cu(001) substrate. The corresponding 
exchange integrals are multiplied by a suitable power of dis- 
tance to highlight their decay exponent. Note that the first 
two exchange integrals in (a) are out of frame. 



acter for the case of partial coverage, namely, decreasing 
number of nearest neighbors which also leads to enhanced 
Fe-moments as compared to a full Fe overlayer. Mathe- 
matically, the exchange integrals are scaled by the same 
factor which also increases the impurity density of states 
(or local moments) as compared to the host one. On 
the other hand, their effective number is proportional to 
x 2 . As a result, the effective interactions (see Eq. ([3])) de- 
crease with decreasing coverage indicating the monotonic 
decrease of the corresponding Curie temperature. 



3 
2.5 

2 
1.5 

1 

0.5 



-0.5 

3 
2.5 

2 
1.5 

1 

0.5 



-0.5 



(a) 



(F ei _ x Co x )ICu(001) 






m 




(b) 






(Fe 50 Co 50 )ICu(001) 
• Fe-Fe 








□ Fe-Co 








O Co-Co 






\ 















0.5 1 1.5 2 2.5 3 3.5 4 
Relative distance D=(d/a) 



FIG. 4: The exchange integrals in (Fea:,Coi_a:) overlayer 
on the fcc-Cu(OOl) substrate as a function of the shell dis- 
tance: (a) The case of pure Fe and Co overlayers and the 
effective VCA exchange integrals ([3]) (Feo.5,Coo.s) overlayer 
(empty squares); and (b) The species-resolved exchange inte- 
grals j9' Q , Q,Q'=Fe,Co for (Fe .5,Co .5) overlayer. 



B. Curie temperature: iron-cobalt overlayer 

We have first investigated the magnetic properties of 
random Fcr^Co^ overlayer on Cu substrate using two 
methods: the RPA and the MC simulation. In MC simu- 
lation we calculated the size-dependent susceptibility and 
then performed extrapolation as described in the subsec- 
tion III B 21 As discussed above, in the RPA we used an 
approximate VCA while MC simulations were done both 
for the VCA as well as for the realistic case with three 
different exchange integrals J®'® (Q,Q'=Fe,Co). The 
RPA- VCA approximation was used successfully for the 
calculation of Tc of bulk Ni-rich transition metal alloys. 21 
In Ref.— , it was demonstrated numerically on a simple 
model that VCA is a good approximation above the per- 
colation limit and for extended exchange integrals. An- 
other reason for validity of the RPA- VCA in this case is 



a similarity of exchange integrals J®'® {Q, Q'=Fe,Co) 
as shown above. On the contrary, for localized exchange 
integrals and a very low concentration of magnetic im- 
purities the VCA fails, like, e.g., in (Ga,Mn)As and 
(Ga,Mn)N diluted semiconductor alloys^ 2 - In the present 
case one should expect that the VCA will be a good ap- 
proximation because of extended character of exchange 
integrals J®'® and their similarity for various atom types 
as indicated above. 

The calculation confirms this prediction as it is ob- 
vious from Fig. [5] in which we compare Tq of fec- 
(Fe,Co)/Cu(001) random overlayer over the whole con- 
centration range. In the limit of Fe- and Co-overlayer on 
fcc-Cu(OOl) the present results agree well with a previ- 
ous study- 2 - and small differences are due to the differences 
in technical details. We have obtained well-pronounced 
maximum in the concentration dependence of Tq . which 
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FIG. 5: The Fe-Fe exchange integrals as a function of the shell 
distance for three different coverages of fcc-Cu(OOl) surface by 
Fe atoms, indicated in the plot. 



reminds a similar maximum in the concentration trend 
of random bulk FeCo alloys in the concentration range 
in which the bcc-phase exists^. The calculated Tc de- 
pends on the values of the anisotropy energy A of con- 
stituent atoms, although such dependence is quite weak 
as demonstrated in Ref. l2|. We remind that this term 
is responsible for the formation of a narrow gap in the 
spin-wave spectra which is an origin of the finite Curie 
temperature for the monolayer case 2 . Due to this term 
the Curie temperature determined in the framework of 
the RPA or Monte Carlo is much lower than its mean- 
field counterpart as it is obvious from Fig. [5] (see also 
Ref. [H) . A reliable choice of this term is a problem. In the 
present case we have chosen the same values as in Ref.! 
(A Fc =0.140 mRy A Co =0.052 mRy) and their average in 
the VCA models in order to have a direct comparison 
with the case of ideal overlayers. A recent studj* 2 ^ indi- 
cates that this value may be significantly smaller. Be- 
cause of this uncertainty, we have tested the robustness 
of the present result with respect to various values of 
this term using the RPA- VCA method. In particular, 
we show in Fig. [6] results of the RPA- VCA for the case 
in which we used the same values of A for both Fe and 
Co, but reduced by an order of magnitude as compared 
to that for Fe-atoms (A=0.015 mRy). Calculated Curie 
temperatures are lowered in agreement with the previous 
studj* 2 -, but the concentration maximum, although less 
pronounced due to the same values of A, is still present. 
It should be noted that a reliable experimental estimate 
of the Curie temperature of monolayers is still a chal- 
lenge. 

We note that similar maxima in the concentration 
trend of Tc exist in random bulk bcc-FeCo^ and fcc-NiFe 
alloys 21 . In the former case the total moment has a weak 
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FIG. 6: (Color online) Concentration dependence of Tc for 
fee (Fe,Co)/Cu(001) random overlayer as a function of Co- 
concentration. Lower panel: We compare the MC simulations 
using random exchange integrals (circles) with the simplified 
MC-VCA (diamonds) and RPA- VCA (triangles and squares). 
In the VCA cases, we employ non-random alloy with effective 
concentration-dependent exchange integrals (see the text). 
Also shown is the RPA- VCA result for reduced anisotropy 
(see the text). Upper panel shows mean field approximation. 



maximum for about 30% of Co2£, in the latter case the 
total moment varies almost linearly with the content 21 
similarly like in the present overlayer case^ This means 
that there is no simple relation between concentration 
trends of total moments and Curie temperatures. The 
Curie temperature is determined by the exchange inte- 
grals. It should be noted that a discussed concentration 
maximum of Tc is seen already in their mean-field val- 
ues, Fig. |51 which, in turn, are directly proportional to 
the sum of averaged exchange integrals, Eq. ([3]). This 
can be considered as a precursor of the concentration 
maximum of Tc . 

We can conclude that a critical comparison of calcu- 
lated Curie temperatures of random magnetic overlay- 
ers on nonmagnetic substrates can be conveniently per- 
formed using a simplified RPA- VCA approach which is 
significantly more numerically efficient than MC simu- 
lations. This is mostly due to the fact that exchange 
integrals for various species are rather similar. Using the 
fast RPA- VCA approach we have verified that due to a 
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slow decay of exchange integrals with distance one should 
carefully check the convergence of results with the num- 
ber of shells included in simulations. 



C. Curie temperature: incomplete coverage by iron 

The next case is the model of partial coverage of the 
fcc-Cu(OOl) substrate by iron atoms. In this case, con- 
trary to the previous model of (Fe,Co)-overlayer, the 
only non-zero exchange integrals are those among Fe 
atoms. The effective exchange integral (VCA) is thus 
JfJ 1 — x 2 J ij e ' e , where x denotes the coverage (x=l corre- 
sponds to ideal Fe-overlayer) . It should be noted that for 
2D-systems the effect of percolation may lead to reduc- 
tion of Curie temperature to zero value. In particular, it 
is known that the site percolation threshold for the square 
lattice with interaction up to next nearest neighbors is 
p c ~ 0.592 and that the percolation threshold may pro- 
gressively decrease with the range of interaction^ pro- 
vided the strength of interaction is independent of a range 
of neighbor bonds. Despite the fact that exchange inte- 
grals are relatively long-ranged, one thus should expect 
that the VCA will represent a substantial approximation, 
at least at low coverages. 

We have estimated the Curie temperature of fcc- 
Fe x /Cu(001) system for coverages x=l, 0.75, 0.5, and 
0.25. The results are summarized in the Table 1. Both 
the RPA-VCA and MC simulations give qualitatively the 
same result, namely, an almost linear decrease of the 
Curie temperature with decreasing coverage. On the 
other hand, the RPA-VCA overestimates the Curie tem- 
perature with decreasing coverage as expected and for 
very low coverage (#=0.25) it gives even the qualitatively 
wrong result as the MC leads to a collapse of the long- 
range magnetic order. One can thus conclude that for 
systems with very different exchange integrals one should 
be careful when using simplified approaches like the RPA- 
VCA. Another example of such system can be, e.g., the 
surface alloy fcc-(Cu,Mn)/Cu(001) with negligible ex- 
change integrals between Cu atoms and also between Cu 
and Mn atoms due to a negligible magnetization of Cu 
atoms, On the other hand, the 2D-generalization of the 
random RPA approach^! which was successfully applied 
to the dilute magnetic semiconductors can be a fast and 
reliable counterpart to the numerically demanding MC 
simulations which, on the other hand, are the most reli- 
able. 



D. Magnon spectra 

Progress in experimental techniques in the recent years 
has made it possible to measure magnon dispersion even 
in ultrathin magnetic overlayers, such as a single Fe 
monolayer on W(110), using spin-polarized electron en- 
ergy loss spectroscopy (SPEELS) measurements^. Al- 
though the present random iron-cobalt overlayer has not 



TABLE I: Curie temperatures of overlayer with incomplete 
coverage by iron. Values of Tc obtained by MC simulation 
and RPA calculation for different concentrations x are shown. 



X 


T C MC (K) 


T C RPA (K) 


1.0 


552 


552 


0.75 


361 


415 


0.5 


140 


269 


0.25 




131 



yet been measured, our calculations hopefully could mo- 
tivate such a study. In principle, there are at least 
three sources of broadening of the magnon dispersion, 
namely, the alloy disorder, the transversal temperature 
fluctuations of the magnetic moments originating from 
the coupled thermal bath, and the longitudinal Stoner 
excitations. At present, the ASD formalism does not 
include Stoner excitations, but for low temperatures 
and smaller wavevectors, the density of magnons dom- 
inates over Stoner excitations. Here, we focus on the 
broadening originating from alloy disorder by calculat- 
ing magnon dispersions for the following cases (i) pure 
Fe, (ii) FesoCoso, and (iii) Fe7sVac25 overlayer on top of 
Cu(001), as displayed in FigJT] The color is a measure 
of magnon damping (the full width at half maximum of 
the magnon spectral function). In all three cases, the 
temperature is fixed at T=5K and a=0.005. We probe 
magnons in the fcc(001) two-dimensional Brillouin zone 
by following the path X ~ f - M - X. The X and M 
points correspond to p(2xl) and c(2x2) antiferromagnetic 
structures, respectively. 

The magnon spectra for the Fe overlayer show very lit- 
tle broadening throughout the Brillouin zone as expected, 
since there is no alloy disorder in this case. Alloying Fe 
with Co changes the situation. Here the spectra are more 
diffuse due to the alloy disorder. It is worth noting that 
the broadening is most pronounced around the M point 
in the Brillouin zone. The higher frequencies at the X 
and M points are consistent with the higher Tc found in 
Feo.sCoo.s overlayer as compared with the pure Fe case. 
The difference between magnetic properties of Fe and Co 
atoms is not very big and there are still distinct magnon 
excitations visible in the whole BZ. However, an extreme 
situation is encountered in the dilute case (Fe75Vac25 al- 
loy). Here the alloy disorder is large causing a very dif- 
fuse magnon spectra. A similar result has been found 
in three-dimensional diluted magnets by random phase 
approximation^. One can say that magnon excitations 
in this case are strongly damped and existing only for 
small wave vectors, i.e., for longwave magnons which are 
less sensitive to local site disorder. In this incomplete 
Fe coverage (x—0.75), we observe a dispersionless exci- 
tation with the energy close to zero. It corresponds to 
localized magnons that appear close and below the per- 
colation threshold. This phenomenon was discussed by 
Chakraborty and Bouzerar ,29] in a three-dimensional 
case, but the same arguments are valid also in two di- 
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mensions. We note that the percolation concentrations 
for fcc-lattice and for the square lattice are around 0.2 
and 0.5, respectively. 



IV. DISCUSSION AND CONCLUSION 

We have presented the first-principles theory of ther- 
modynamical properties of random magnetic overlayers 
on non-magnetic metallic substrates and applied it to two 
systems, namely, to a random iron-cobalt monolayer on 
the fcc-Cu(OOl) substrate and to the case of partial cov- 
erage of fcc-(OOl) face of copper by iron atoms. Atomistic 
spin dynamics simulations were used to predict magnon 
dispersions in random overlayers, where broadening from 
alloy disorder was quantified. 

The main conclusions from our study are: (i) The ex- 
change integrals of magnetic overlayers decay more slowly 
as compared to the bulk cases and large number of shells 
has to be included in the statistical treatment to obtain 
well converged result; (ii) Contrary to the bulk case, we 
have found that the most efficient and reliable approach 
to estimate the Curie temperature of magnetic overlay- 
ers is the estimate of the susceptibility combined with 
the scaling relation rather than the cumulant method; 
(iii) The simplified, but numerically efficient RPA-VCA 
approach was tested and limits of its applicability were 
established. It was found that the RPA-VCA is applica- 
ble for systems with similar exchange integrals (like, e.g., 
FeCo-overlayer) , but it fails for systems with very dif- 
ferent values of exchange integrals (like, e.g., incomplete 
Fe-overlayer) ; (iv) Calculations indicate the presence of a 
maximum in the concentration dependence of the Curie 
temperature similar to that observed in bulk bcc-FeCo 
alloys. 
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